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Taking advantage of an exact mapping between a relativistic integrable model and the bieb-biniger 
model we present a novel method to compute expectation values in the bieb-biniger Bose gas both 
at zero and finite temperature. These quantities, relevant in the physics of one-dimensional ultracold 
Bose gases, are expressed by a series that has a remarkable behavior of convergence. Among other 
results, we show the computation of the three-body expectation value at finite temperature, a 
quantity that rules the recombination rate of the Bose gas. 
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Correlation functions are key quantities in quantum 
interacting systems: not only they fully encode the dy- 
namics but also are directly related to various susceptibil- 
ities and response functions. For these reasons, there has 
always been an intense search to find the most efficient 
ways to compute them. The task is notoriously difficult, 
even if the system is integrable. A typical but signifi- 
cant example is provided by the bieb-Liniger (LL) model 
PQ that describes the low-temperature properties of one- 
dimensional interacting Bose gases: although it can be 
solved through Bethe Ansatz equations, the explicit com- 
putation of its correlation functions is a long-standing 
problem [2j [3] . The interest for the computation of cor- 
relation functions in the bb model is obviously not only 
theoretical. In a series of recent experimental achieve- 
ments strongly interacting ultracold bosons have been 
confined within wave-guides by nearly one-dimensional 
potentials that tightly trap the particle motion in the 
two transverse directions while leaving it free in the third 
axial direction [IHE|: the coupling of these bosons with 
the external world can be made so weak that their be- 
haviour is very well described by the bb model. Through 
interference (or eventually in situ) experiments, several 
quantities can be detected both at zero and finite temper- 
ature: the time duration of experiments depends on the 
three-body recombination rate, which is proportional to 
local three-body expectation values [7] . Many important 
general questions of quantum many-body physics can be 
studied in such a highly controllable set-up: dynamical 
properties concerning the absence (or not) of thermaliza- 
tion jH [5], for instance, or the behaviour of integrable 
quantum systems when small non-integrable perturba- 
tions (e.g., three-body interactions and/or a weak exter- 
nal trapping potential) are switched on [S]. 

Over the years several theoretical quantities of the bb 

model have been computed by means of different tech- 
niques |f 0H19j . In this paper we present a compact and 
general way to determine the expectation values of its 
local operators. The method takes advantage of an exact 
mapping between a relativistic integrable massive model 
- the sinh-Gordon (ShG) - and the bb model: in a proper 
non-relativistic limit of the ShG model, both its S'-matrix 
and bagrangian coincide with those of the bb model. 



Since the S'-matrix of an integrable relativistic model 
fixes the exact matrix elements of all operators of the 
theory (and for the ShG model these matrix elements are 
all known [201 [2T]), the correspondence between the two 
models opens the way to computing the corresponding 
quantities of the bb model in a very direct way. As shown 
below, this method provides a remarkable simplification 
of the problem. Its implementation actually requires to 
take into account an additional aspect of the problem: 
while in the ShG model the correlation functions refer 
to the vacuum (i.e. the state without any particles), in 
the bb model they relate instead to its ground state at 
a finite density. This aspect, however, can be success- 
fully overcome by the Thermodynamical Bethe Ansatz 
(TBA) formalism developed in [52], which has the addi- 
tional convenience of being applicable equally well both 
at zero and finite temperature. In this way we are able 
to compute not only the zero temperature expectation 
values but also their finite temperature expressions. 

The LL-ShG mapping. The bb Hamiltonian for N in- 
teracting bosons of mass m in one dimension is 

h 2 N Fi 2 

i— 1 1 i<j 

The corresponding non-relativistic field theory descrip- 
tion is the quantum non-linear Schrodinger model [2], 
which employs the complex field ip and the bagrangian 

The effective coupling constant of the bb model is given 
by the dimensionless parameter 7 = 2m\/h 2 n, where 
A > is the coupling entering the Hamiltonian while 
n = N/L is the density of the gas (L is the length of 
the system) . Temperatures are usually expressed in units 
of the temperature Tr> = h 2 n 2 /2mfcs of the quantum 
degeneracy, r = T/Trj. The two-body elastic S-matrix of 
the bb model is [3 [TO] 

Sll (P, A) = , (3) 
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where p is the momentum difference of the two particles. 

Consider now the ShG model in (1+1) dimensions, i.e. 
the integrable and relativistic invariant field theory de- 
fined by the Lagrangian 
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where <fi = <j)(x, t) is a real scalar field, fi is a mass 
scale and c is the speed of light. The parameter [i is re- 
lated to the physical (renormalized) mass m by /i 2 = 
-Kam 2 c 2 /h 2 sin(7ra), where a = he g 2 / (8tt + he g 2 ) [20] . 
The energy E and the momentum P of a particle can be 
written as E = mc 2 cosh#, P = mcsmh.9, where 9 is the 
rapidity. Since the ShG dynamics is ruled by an infinite 
number of conservation laws, all its scattering processes 
are purely elastic and can be factorized in terms of the 
two-body S-matrices [20] 



<SshG(0,a) 
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sinh 9 + i sin(a7r) 
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where 9 is the rapidity difference of the two particles. It 
is now easy to see that taking simultaneously the non- 
relativistic and weak-coupling limits of the ShG model 
such that 
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its S'-matrix ([5| becomes identical to the S'-matrix ([3| 
of the LL model. Notice that the coupling A does not 
need to be small, i.e. with this mapping we can study the 
LL model at arbitrarily large values of the dimensionless 
coupling 7. 

The mapping between the two models goes beyond the 
identity of their S-matrix: it actually extends both to 
their Lagrangians and TBA equations. Details will be 
given elsewhere , but it is simple to follow the main steps 
of the procedure. According to [23], the non-relativistic 
limit of a field theory consists of expressing the real scalar 
field in the form 



{x,t) = 




+ ^ t (.T,t)e +4 ^ t ) , 



and, when the limit c — > 00 of the Lagrangian is taken, 
of omitting all its oscillating terms . The commutation 
relation [<fi(x, t), Tl(x', t)} = ihd(x — x') implies for the 
non-relativistic operators [t/j(x,t),ip^(x',t)] = S(x — x') . 
Furthermore, when the limit g — > of eqn ^ is con- 
sidered, the ijftij) terms coming from the potential and 
kinetic parts cancel each other, while all higher terms 
of the series expansion of the potential, but the quartic 
one, vanish. Hence, the ShG Lagrangian Q reduces to 
the non-linear Schrodinger Lagrangian Notice that 
the mapping based on the limit ^ applies to any oper- 
ator of the theory. 

In the same way one can also show that the TBA equa- 
tions of the ShG model (given for instance in [53] ) reduces 



to the ones of the LL model, written down in [ID] . In the 
LL model at a finite T the TBA equation for the pseudo- 
energy e(T, fi) consists of the non-linear integral equation 



p 2 /2m — /1 
k B T 



if o log (l + e 
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where fj, is the chemical potential associated to the finite 

log Sll (p) is the deriva- 



density n of the gas, ip(p) — — t,-^ 

tive of the phase shift and ipof = <p(p—p')f(p') 
The solution of this integral equation leads to the free en- 
ergy and to all other thermodynamical data of the model. 

Expectation values. At equilibrium the expectation 
value of an operator O = 0(x) at temperature T and 
at finite density is given by 
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In a relativistic integrable model the above quantity can 
be neatly expressed as [52] 
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where f(9i) = 1/(1 + e € ^>) and 9 = x ,...,6 n 
(9 = 9 n ,...,9i) denote the asymptotic states enter- 
ing the traces in This formula employs both the 
pseudo-energy e(9) and the connected diagonal form 
factor of the operator O, defined as (9 \0\ 9 ) conn = 



i^P(lim ?Ji _,o(0|e>| 9' , '9 -ivr + iV)) where V = Vn, ■ ■ ■ ,Vl 
and FP in front of the expression means taking its finite 
part, i.e. omitting all the terms of the form r]i/r)j and 
1/rif where p is a positive integer. 

In order to compute the expectation values of the LL 
model by applying eqn ^ we need: (a) to solve the inte- 
gral equation Q for e(0); (b) to identify the relevant form 
factors of the ShG model; (c) to take the non-relativistic 
limit of both the form factors and eqn Given for 
granted the straightforward (numerical) solution of eqn 
([7]), let us focus our attention on the last two points. The 
generic m-particle form factor of a local operator O in 
the ShG model can be written as [20] [H] 
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where 2^ = e 6i and Q® n are the symmetric polynomials 
in the x's that fully characterize the operator O. The 
explicit expression of F m i n (9) is given in [5D] but the only 
thing needed here is its functional equation 

sinh 9 

F min (™ + 0)F min (9) = - 



sinh 9 + i sin(7ra) 



We are interested in the symmetric polynomials Qm' of 
the exponentials O q = e q9 ^ since, using their Taylor ex- 
pansion, we can extract the form factors of all normal 
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FIG. 1: gi and g2 at T = using form factors up to n = , 
6 and 8 particles, respectively with green dot-dashed, blue 
dashed and red dotted lines. The exact values are given by 
the solid lines whereas the dot-dashed line below, indicated 
by the arrow, corresponds to the strong coupling expansion 
dl5». 



; : . Their expression is 
4sin(7ra) 



ordered operators : 

Q# = [q] 



N 



detM m (<j), (11) 



where M m (q) is an (m — 1) x (m — 1) matrix with el- 
ements [M m (q)] itj = viT-jii ~ 3 + Q]- Above, [x] = 
sin(a;7ra)/ sin(7ra) while ci™ 1 ' (a = 0,1,..., m) are the 
elementary symmetric polynomials in m variables. 

On the basis of the results given above, we are now in 
the position to compute the local /c-particle correlation 
function g k of the LL model denned by 



(^V) =n k g k ( 1 ,r) 



(12) 



where k is an integer (k = 1, 2, 3, . . . ). The g k 's are func- 
tions of the dimensionless LL coupling 7 and of the re- 
duced temperature r. The relation between g k in the LL 
model and the corresponding quantity in the ShG model 
in the limit ^ is given by 
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Using eqn (JoJ) and the connected form factors of the cor- 
responding operator we arrive at the expression 
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are the double limit Q of the connected form factors. As 
shown below, the series ( 13 ) are nicely saturated by the 



FIG. 2: gi at r = 1, 10 using form factors up to n — L, 6 and 
8 particles with green dot-dashed, blue dashed and red dotted 
lines, respectively. The solid lines show the exact result. 



first few terms for sufficiently large values of 7 (7 = 
is a singular point of the model [T], therefore one can- 
not expect a priori any fast convergence nearby). A first 
check of eqn (13) is provided by the case k = 1: using 



( 13 1 (with a chemical potential fi that ensures the finite 



density n) and summing up the series, one easily checks 
that (tp^) — n an d 9i = 1, as it should be for trans- 
lational invariance. As shown in Fig. [T] the exact value 
gi = 1 (solid line) is rapidly approached by just the first 
terms of (13 1: the convergence of the series is always re- 



markably fast for all 7 > 1.5, where the exact value is 
obtained within a 5% accuracy just using its first four 



terms. As a second check of ( 13 1 let us show how we can 



easily recover the leading order of the strong coupling 
(i.e. large 7) expansion of all g k - since this always comes 



from the first non-zero integral in the series (13), we get 
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where I k = J_ 1 dfci ■■•J_ 1 dk k JJ^jih — kj) 2 . This result 
coincides with the one obtained in [14]. 

The quantity g 2 can be exactly determined via the 
Hellmann-Feynman theorem |15j and its plot at T = is 
shown in Fig.[l]together with our determination from eqn 
( 13 ). As before, also in this case there is a fast convergent 



behaviour of the series. The strong coupling regime of 52 
can be computed by expanding in powers of j^ 1 all the 



terms in eqn ( 13 ) and for T = we get 
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in agreement with the result of the Hellmann-Feynman 
theorem [THUS]. Expression (15) is also plotted in Fig.[l] 
in order to show that the determination of <?2 (at finite 
7) obtained from the first terms of eqn (13) is closer to 



the exact result, because any of them contains infinitely 
many powers of 7. At finite temperatures the convergence 
of the series is also pretty good and the results are shown 
in Fig. [2} 

As a final example, let us discuss g^, a quantity known 
exactly at T = [TH], but only approximately at T > 
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FIG. 3: 33 vs 7 at : (a) r = 0, (b) r = 1 and (c) r = 10. In (d) 
we plot 33 vs r at 7 = 7. In all figures the blue dashed and the 
red dotted lines refer to n = 6 and 8 particles, respectively. 
The solid line in (a) is the exact value of 33 at r = 0. 



[13] . From (13) its strong coupling limit at T = is 



93 



16 7T 6 
I576 



16 

1 



+ 0(7~ 8 ) 



(16) 



where we report both the leading and sub-leading terms 
in 7 _1 of this expression. The plot of 33 at t = using 
form factors up to n — 6 and 8 particles (i.e. one or two 
terms of the series ( 13 )) is in Fig. [3^a) and, as in previous 
examples, it shows a nice convergent pattern to the exact 
value found in [19] . Figs. [3]^b,c) show 33 as a function of 
7 at fixed temperature r, while Fig. [3jci) shows instead 
33 as a function of r at a fixed value of 7. 

Conclusion. We have shown that the equilibrium ex- 
pectation values for one-dimensional interacting Bose 
gases can be efficiently computed by using the non- 
relativistic limit of an integrable relativistic field theory, 
the sinh-Gordon model. There is a significant advantage 
in using this method instead of employing directly the 
non-relativistic Lagrangian {2). The reason is that a rel- 
ativistic field theory presents a larger number of con- 
straints (crossing invariance, for instance) which permit 
to pin down exactly and efficiently the matrix elements 
of all operators: once these quantities are known, it is 
then easy to take their non-relativistic limit. As shown 
above, this proves to be a notable simplification in the 
computation of the correlators of the LL model. 

The method works equally well at T = and T ^ 
where the series expansion presents a remarkable conver- 
gence behaviour for finite values of 7. There is no ob- 
struction, in principle, to compute higher form factors 
and further improve the result. Strong coupling expan- 
sions in 7 _1 can be easily derived as well but at finite 
7 the form factor expansion, containing infinitely many 



powers of 7, is more accurate for the determination of gk, 
as we showed comparing it with exact results. As a sig- 
nificant application of the method we have determined 33 
at finite temperature (a term which is proportional to the 
recombination rate of the gas) . This quantity, as well as 
the higher g^, may provide important information once 
the integrability of the model is broken. In the future it 
would be also interesting to apply this method both to 
two-point correlation functions and to other models. 
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